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Abstract 

The theory of turbulent photon filamentation in lasers with high Fresnel num- 
bers is presented. A survey of experimental observations of turbulent filamentation is 
given. Theoretical description is based on the method of eliminating field variables, 
which yields the pseudospin laser equations. These are treated by the scale separation 
approach, including the randomization of local fields and the method of stochastic av- 
eraging. The initial, as well as the transient and final stages of radiation dynamics are 
carefully analysed. The characteristics of photon filaments are obtained by involving 
the probabilistic approach to pattern selection. 
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1 Introduction 



In nonlinear media interacting with electromagnetic fields there appear different spatiotem- 
poral structures that are analogous to the structures arising in many other complex nonequi- 
librium systems [1-5]. The most known among such electromagnetic structures are the opti- 
cal filaments which can be formed in passive nonlinear matter [3,4,6] and in active laser media 
[1,4,5]. Examples of filaments occurring in other nonequilibrium systems can be found in 
the reviews [2,5]. In addition, it is possible to mention the generation of filaments in hadron 
plasma [7] and in quark-gluon plasma formed under heavy nucleus collisions [8,9]. 

The present survey will be focused on the optical filaments developing in lasers. The 
behaviour and characteristics of these filaments essentially depend on the value of the Fresnel 
number F = R 2 /XL, in which R and L are the internal radius, aperture radius, and effective 
length, respectively, of a cylindrical laser, and A is the optical wavelength. There are two 
types of optical filaments, regular and turbulent, corresponding to either low or high Fresnel 
numbers. The difference between these two filament types, their main properties, and the 
related experiments are discussed in Sec. 2. 

Our concern here is the turbulent filaments arising in lasers with high Fresnel numbers. 
The transition from a regular filamentary structure to the turbulent filament behaviour is 
somewhat analogous to a crossover phase transition in statistical systems [9-11], or a more 
close analogy is the transition from the laminar to turbulent motion in liquids [1,2,12]. The 
description of the turbulent motion is notorious to be rather difficult and is usually done by 
invoking a kind of averaging procedure. Optical turbulent structures also require the usage 
of averaging techniques. 

First, the turbulent filamentation in laser media was described on the basis of stationary 
models [13-16] invoking the notion of an effective time-averaged energy. A more elaborated 
approach, based on realistic evolution equations, was developed in Refs. [5,17-20]. It turned 
out that for this purpose it is convenient to work not with the standard Maxwell-Bloch 
equations but with the evolution equations for the pseudospin operators, resulting after the 
elimination of the field variables. Since this type of equations is less known for collective 
electromagnetic phenomena, these equations are derived in Sec. 3, and their specification 
for treating the turbulent filamentary structure is given in Sec. 4. 

Turbulent filamentation in lasers is a self-organized process due to the photon exchange 
through the common radiation field. Coherent radiation in filaments develops in a self- 
organized way, even when there has been no coherence at the initial time. The triggering 
mechanism for the origination of the radiation coherence is the existence of transition dipolar 
waves caused by the dipolar part of effective atomic interactions. This is emphasized in Sec. 
5. 

The overall dynamics of radiating filaments is analysed in Sec. 6, starting from the 
quantum stage, through the transient flashing regime, to the stationary state. Generally, 
the filaments of different radii can emerge, whose classification is given by means of the 
probabilistic approach to pattern selection [18,21,22]. This allows us to define the typical 
filament radius and the number of filaments, which is done in Sec. 7. 

The present work does not merely summarize and review the previous theoretical results 
on the turbulent filamentation in laser media but also contains several improvements of the 
theory, making the latter better grounded both mathematically and physically. 
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2 Experiments on Filamentation in Lasers 



First of all, let us recall that there are two types of optical filaments, regular and turbulent, 
which is related to the value of the laser Fresnel number. The latter plays for optical 
systems the same role as the Reynolds number for moving fluids. When increasing the 
Reynolds number, a laminar fluid transforms to turbulent. In the similar manner, increasing 
the Fresnel number makes a regular filamentary structure turbulent. Optical turbulence 
implies, by analogy with the fluid turbulence, that the spatiotemporal dynamics is chaotic. 
This means that the radiating filaments are randomly distributed in space and are not 
correlated with each other. 

Optical filaments are observed in the near-field cross-section of lasers. The typical pic- 
ture, when varying the Fresnel number is as follows. At very small Fresnel numbers F <^1, 
there exists the sole transverse central mode uniformly filling the laser medium. When the 
Fresnel number is around F ~ 1, the laser cavity can house several transverse modes seen 
as a regular arrangement of bright spots in the transverse cross-section. Each mode corre- 
sponds to a filament extended through the cylindrical volume. This filamentary structure 
is regular in space, forming ordered geometric arrays, such as polygons. The transverse 
structure is imposed by the cavity geometry, being prescribed by the empty-cavity Gauss- 
Laguerre modes. Such regular structures are well understood theoretically, their description 
being based on the field expansion over the modal Gauss-Laguerre functions related to the 
cylindrical geometry [1]. For Fresnel numbers up to F m 5, the number of bright filaments 
follows the F 2 law as F increases. The regular filamentary structures have been observed 
in several lasers, such as CO2 and Na2 lasers [1]. Similar structures also appear in many 
passive nonlinear media, e.g. in Kerr medium [3,4] and in active nonlinear media, as the 
photorefractive Bii 2 Si0 2 crystal pumped by a laser [3,4]. 

As soon as the Fresnel number reaches F pa 10, there occurs a qualitative change in 
the features of the filamentary structure: The regular filaments become turbulent. This 
transition goes gradually, as a crossover, with the intermittent behaviour in the region 5 < 
F < 15. The character of this change is again common for lasers [23,24] and for active 
nonlinear media [3,4]. 

For Fresnel numbers F > 15, the arising filamentary structures are principally different 
from those existing at low Fresnel numbers. The spatial structures now have no relation to 
the empty-cavity modes. The modal expansion is no longer relevant and the boundary condi- 
tions have no importance. The laser medium houses a large number of parallel independent 
filaments exhibiting themselves as a set of bright spots randomly distributed in the trans- 
verse cross-section. The number of these random filaments is proportional to F, contrary to 
the case of low Fresnel numbers, when the number of filaments is proportional to F 2 . The 
chaotic filaments, being randomly distributed is space, are not correlated with each other. 
Such a spatio-temporal chaotic behaviour is characteristic of hydrodynamic turbulence, be- 
cause of which the similar phenomenon in optics is commonly called the optical turbulence. 
In contrast to the regime of low F, where the regularity of spatial structures is prescribed 
by the cavity geometry and boundary conditions imposing their symmetry constraints, the 
turbulent optical filamentation is strictly self-organized, with its organization emerging from 
intrinsic properties of the medium. Since the optical turbulence is accompanied by the for- 
mation of bright filaments with a high density of photons, this phenomenon can be named 
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[17] the turbulent photon filamentation. This phenomenon is common for lasers as well as 
for photorefractive crystals [3-5]. 

The first observations of the turbulent filamentary structures in lasers, to my knowledge, 
were accomplished in the series of experiments [25-29] with the resonatorless superlumines- 
cent lasers on the vapours of Ne, Tl, Pb, N 2 , and N^. In these experiments, the typical 
characteristics were as follows: A pa 5 x cm, R pa 0.1 — 0.3 cm, L pa 20 — 50 cm, and 
F pa 10 — 100. The number of filaments was Nj ~ 10 2 — 10 3 , with the typical radius 77 pa 0.01 
cm. 

Then the filamentary structures in large-aperture optical devices have been observed 
in several lasers, as reviewed in [23,24], and in photorefractive crystals [3,4]. Numerical 
simulations have been accomplished [30]. Experimental works mainly dealt with the CO2 
lasers [23,24,31,32], dye lasers [33], and semiconductor lasers [34,35]. 

The turbulent nature of filamentation occurring in high Fresnel number lasers was care- 
fully studied in a series of nice experiments [36-42] with C0 2 lasers and dye lasers. Irreg- 
ular temporal behaviour was observed in local field measurements. It was found that the 
transverse correlation length was rather short. Randomly distributed transverse patterns 
generated in short times were observed, being shot-to-shot nonreproducible. For intermedi- 
ate Fresnel numbers F ~ 10, instantaneous transverse structures were randomly distributed 
in space, but after being temporally averaged, they displayed a kind of regularity related to 
the geometrical boundary conditions. This type of combination of irregular instantaneous 
patterns with the averaged or stationary pattern, showing the remnant ordering, is under- 
standable for the intermediate regime in the crossover region 5 < F < 15. Fully developed 
optical turbulence is reached as the Fresnel number increases up to F ~ 100. 

The typical laser parameters are as follows [36-42]. The pulsed C0 2 laser, with the 
wavelength A = 1.06 x 10~ 3 cm and frequency 00 = 1.78 x 10 14 s -1 , emits the pulses of 
t p pa 0.7 x 10 -7 s or 10~ 6 s. The aperture radius R pa 1 cm, laser length L = 100 cm. The 
inversion and polarization decay rates are 71 = 10 7 s _1 and 72 = 3 x 10 9 s -1 . The CO2 
density is p = 2 x 10 18 cm~ 3 . The Fresnel number is F pa 10. The characteristic filament 
radius is 77 pa 0.1 cm. 

The pulsed dye laser, with the wavelength A = 0.6 x 10~ 4 cm and frequency uo = 3.14x 10 15 
s _1 , produces pulses of t p pa 0.5 x 10~ 6 s. The decay rates are 71 = 4 x 10 8 s" 1 and 72 = 10 12 
s -1 . The cavity length is L pa 20 cm. By varying the aperture radius between 0.3 cm and 
0.8 cm, the Fresnel number can be changed by an order, between F = 15 and F — 110. The 
typical filaments radius is 77 pa 0.01 cm. 

In a broad-aperture pulsed dye laser, two different spatial scales were noticed [33], one 
order of magnitude apart. The appearance of the second, larger, scale could be due to a 
nonvanishing interaction between filaments. 

3 Pseudospin Laser Equations 

Spatiotemporal laser dynamics is usually treated in the frame of the Maxwell-Bloch equations 
[43]. One commonly employs these equations for defining the points of stability of uniform 
solutions, which means the appearance of nonuniform structures. The linear stability analysis 
of the Maxwell-Bloch equations is the maximum one can do in the case of well developed 
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optical turbulence. 

Another approach is based on the elimination of field variables and dealing solely with 
the equations for spin operators [5,44]. This approach possesses the following advantages: (i) 
It is microscopic, which allows us to better understand the underlying physics of collective 
effects, (ii) It is quantum, which allows for the description, when coherence has not yet been 
developed, (iii) It provides us the possibility not solely for finding instability points but for 
describing the whole dynamics under the condition of strong turbulence. Since this approach 
is not widely known, the related evolution equations are derived below. 

Aiming at pursuing a microscopic picture, let us start with a realistic Hamiltonian 

H = H a + H f + H af + H mf (1) 

representing a system of resonant two-level atoms plus electrodynamic field. The atomic 
Hamiltonian is 

N /l \ 

ifa = I> (2 + 5 *) > ( 2 ) 

where N is the number of atoms; ujq, carrying transition frequency; Sf, spin operator of an 
i-th atom. Wishing to be more rigorous in terminology, we should call the operators Sf, 
appearing here and in what follows, the pseudospin operators, since they correspond not 
to actual spins but to the population difference and dipole transition operators. But for 
brevity, one often calls them just spin operators, which should not bring confusion. The field 
Hamiltonian is 

^ = ^/( E2+H V r < 3 > 

with electric field E and magnetic field H = V x A. The vector potential A is assumed to 
satisfy the Coulomb calibration 

V • A = . (4) 
The atom-field interaction is described by the Hamiltonian 

N /I \ 
#«/ = -£ -Ji-Ai + Pi-Eoi , (5) 



in which the short-hand notation is used for the vector potential Aj = A(rj,t) and the 
external electric field E j = E (i"j,t), and where the transition current is 

J, = iu (dS+ - d*S~) (6) 

and the transition polarization is 

Pi = dSf + d*Sr, (7) 

with the atomic transition dipole d and the ladder operators Sf = S*±iSf. Except resonant 
atoms, the laser cavity can contain some additional filling matter [45] interacting with the 
electromagnetic field through the Hamiltonian 

H m f = ~~ c J}mat(r, t) • A(r, t) dv , (8) 
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where j mai is the local density of current in the filling matter. Here and in what follows, we 
set h = 1. 

The field operators satisfy the equal-time commutation relations 



E a (r,t), AV,*)] = 47ncM r - r ') > 
£ Q (r,t), #V,*)] = -4«c^e a/37 A 5(r - r') , 



(9) 



in which e Q/ g 7 is the unitary antisymmetric tensor [46] and the so-called transverse 5-function 
is 

k a kP\ ik . r dk 



<W(r) = / [S a 3- 



M(') + ^ / V (2^ = I ^ 5(r) - i D ^ (r) ' 
where in the dipolar tensor 



k 2 



(27T)* 



3^^ 



DaetT) 

one has n = r/r = and r = |r|. 

The Heisenberg equations of motion for the field operators yield 

- — - V x H — ' - — -E 

c dt c ^ ' c <9t 



(10) 



;i2) 



from where, with the Coulomb calibration (4), one has the equation for the vector potential 

(13) 



; 2 i d 2 \ A 4vr . 

V -~ 2 W 2 A = "T J ' 



with the density of current 



/3 



- r,)jf (t) + / ^(r - r')j'L t (r',t) dr' 
The known solution to Eq. (13) is the sum 
A(r,*)=A TOC (r,i) + i 



(14) 



r — r 



dr' 



r — r' 



(15) 



of the vacuum vector potential and the retarded potential. 

The Heisenberg equations of motion for the spin operators, satisfying the commutation 
relations 

S~] = 25^ , \Sf, Sf] = ±^-5f 



lead the equations 



dSj 
dt 



= -iu Sr + 25? (fc d ■ A t - id ■ E 0l ) , 
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dS[ 
dt 



-S+ (k d • Aj - id • E 0l ) - S~ (k d* ■ Ai + td* • Eoi) . 



These are complimented by the retardation condition 

S?(t) = {t < 0) . 



(16) 



(17) 



In order that the notion of resonant atoms would have sense, one requires, as usual, 
that the atom-field interactions are small compared to the atomic transition energy, so that 
|d • E| <C lu . Because of this, the retardation can be taken into account in the Born 
approximation 



ST It 



S-(t)Q(ct-r) e 



ikor 



s*lt 



SUt)e(ct-r), 



(18) 



where k = u /c and Q(t) is the unit-step function. 

The idea of eliminating the field operators is based on the usage of the pseudospin 
equations (16), with the substituted there the vector potential (15). In this way, one meets 
the terms corresponding to the atomic self-action, which in the present approach can be 
treated as follows. The vector potential generated by a single atom is 



(19) 



with the current 



J t- 



iu 



dS + {t) e~ lk ° r - d* S~(t) e tkor 6(ct - r) , 



where S a (t) = ^"(O, t). At small distance, such that k r <^ 1, one may write e lk ° r ~ 1 +ik r. 
Substituting the transverse 5-function (10) into the vector potential (19), we keep in mind 
that averaging the dipolar tensor (11) over spherical angles gives 



/ 



D a/3 (r) dn(r) = . 

Then, for k r 1, the vector potential (19) becomes 

2k 



A s (r,t) ~ -k 2 [dS+(t) +d*S~(t) 



+ i 



3r 



dS + {t) -d*S~(t) 



(20) 



To avoid the divergence in the term 1/r, let us average it between the electron wavelength 
A e = 27ih/mc, with m being the electron mass, and the radiation wavelength Ao = 2n/ko. 
Taking into account that A e Aq, we have 



1 



A d r 

e r 



Ao — A e 

Then for the self-acting vector potential, we get 



ko , ( m c 
— In 

2vr 



hcoo 



A s (0,t) = I kl [dS+(t) + d*S-(t)] + ^ In [dS+(t) - d*S~(t) 



2\ 



(21) 
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Substituting this into Eqs. (16) for the case of a single atom, we employ the properties of 
operators of spin 1/2, 

s-s- = s + s + = o , s z s z = \ , s~s z = \s- , s z s- = - \ s~ , 
S + S z = -- s + , s z s + = - s + , s-s + = - - s z , s+s- = - + s z . 

2 2 2 2 

Then we come to the equations for a single atom 

^ = -i(wb -5 L - i lQ )S~ + ^ (7o + iS L )S + , ^ = -2 7o (1 + 5*) , (22) 

in which the notation for the natural width 

7o = ^ |d| 2 fc 3 (23) 

and the Lamb shift 

(24) 

are introduced. The solutions to Eqs. (22), keeping in mind that 70 <C c<Jo and 5l <C are 

S - ^) = S~(0) exp {-i{uj - 5 L )t - 7 t} , 



s z (t) = - 1 - + 



l + s z (o) 



exp(-27 t) . 



Thus, the existence of self-action leads to the appearance of attenuation in the dynamics 
of the spin operators and to the Lamb frequency shift. The latter can always be included 
in the definition of the transition frequency c^o- Taking into consideration the attenuation, 
one usually generalizes the equations of motion by including 72 , instead of 70, for S~ and 
inserting 71, instead of 270, for S z . 

For a system of iV radiating atoms, the vector potential (15) can be presented as a sum 

Here A vac is caused by vacuum fluctuations. The vector potential 

is due to the spherical part of the potential (15), produced by radiating atoms, which in 
addition produce the dipolar part 

W) = -EE/%^4-^)*<. m 

Finally, the action of matter, filling the cavity, creates the vector potential 

A a mat (r, t) = £ / 5a ^'~*P jLt ( r", t - J^l) dr'dr" . (28) 
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Let us combine the vacuum, dipole, and matter vector potentials into the sum 



£(r, t) = 2k d ■ (A vac + A dip + A mat ) , (29) 

which describes local field fluctuations. Owing to the local nature of the fluctuating field 

(29) , it can be treated as a random variable. Contrary to this, the radiation potential (26) 
is of long-range nature and can be responsible for collective effects. The existence of two 
types of variables, acting on different spatial scales, makes it possible to employ the scale 
separation approach [5,47-50]. Then A rad and £ are considered as different operators. Since 
A rad is expressed through the spin operators, we may treat the set S = {Sj\ j — 1,2, ... , N} 
of these operators as approximately commuting with £. Thus, the total set of operators 
consists of two types of the operators, S and £. For any operator function F = F(S,£), we 
may introduce two kinds of averages. One is the average over the spin variables, 

<F> = Tr s pF(S,0, (30) 
with the trace over spins, and p being a statistical operator. Another average is defined as 

«F» = Tr^F(S,0, (31) 

with the trace over the stochastic field (29). 

For the spin averaging (30), we may use the decoupling 

<S?Sf> = <S?><Sf> (i^j), (32) 

keeping in mind the long-range nature of A rad . This reminds us the mean-field approxi- 
mation, which is valid for long-range forces [51]. However, there is a principal difference 
between the mean-field approximation and the decoupling (32). The latter involves only the 
spin degrees of freedom, not touching the stochastic variables £, which are responsible for 
quantum effects. Employing a seemingly semiclassical form (32), at the same time preserv- 
ing quantum features, associated with stochastic field £, is close to the method of stochastic 
quantization used in quantum field theory [52]. Therefore the decoupling (32) can be called 
the stochastic mean-field approximation or the quantum mean-field approximation. 

Let us now average the operator equations (16) over the spin variables, according to Eq. 

(30) , with employing the following notation. The transition function 

u l {t) = u{v l ,t) = 2<Sr{t)> (33) 

describes the local dipole transitions. The local characteristic of coherence is the coherence 
intensity 

2 ™° r i 

Wi (t) = wfat) = - E < S?(t)ST(t) > + < S+(t)Sr(t) >] , (34) 
n ° <j> 

in which the sum is over the nearest neighbours and no is the number of the latter. Finally, 
the local population difference is 

s i (t) = s(r i ,t) = 2<S?(t)> . (35) 
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It is also convenient to pass to the continuous spatial representation, replacing the sums 
by the integrals as 

. , £)■ (») 



A? 



Let us introduce the notation for the effective field potential acting on atoms, 

/(r,t) = / (r,t) + / rad (r,t)+£(r,t), (37) 

which consists of the part 

/ (r,t) = -2id.E (r,i), (38) 
due to an external electric field E , of the term 

f rad (r, t) = 2k <d- A rad {r, t) > , (39) 

caused by the radiating atoms, and of the fluctuating random field (29). The radiation 
potential (39), taking account of Eq. (26), can be written as 



frad(r,t) = -ry p J 
with the transfer kernel 



G(r - r', t)u(r', t) - j^G*(r- r', t)u*(r\ t) 



e i q^ 
k n r 



dr' , 



In this way, from Eqs. (16), we obtain for the local functions (33) to (35) the evolution 
equations 

= -(iLU + -/ 2 )u + fs , — = -272W + (u*f + f*u)s , 



dt 



dt 

% = ~\ («V + /*«)- 7i(« - , 



(40) 



where 71 and 72 are the longitudinal and transverse attenuation rates and £ is the stationary 
pumping parameter. These stochastic differential equations are the basic equations to be 
used in what follows for describing the spatio-temporal evolution in lasers. 



4 Turbulent Photon Filamentation 

The external field E is here the seed field 

E (r, t) = - Eie i(fc * -W * ) + - EJe - *^^ , (41) 

selecting a longitudinal mode of frequency u> = kc, but imposing no constraints on possible 
transverse modes. The propagation of the field (41) is along the axis z, which is the axis of 
a cylindrical laser cavity. The frequency u has to be in resonance with the atomic transition 
frequency u>o, so that the detuning be small, 

1^ < 1 (A = u - u Q ) . (42) 
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The wavelength A = 2-kc/uj is usually much smaller than the effective laser radius R and 
length L, 

A « 1 , A « i . (43) 
The seed field (41) is very weak, such that 

< 1 (i/i = d • Ei) . (44) 

Keeping in mind the possibility of arising filamentary structures, we should look for the 
solution of Eqs. (40) in the form of the modal superposition 

N f N f 



u(r, t) = J2 Mr±, t)e ikz , w(r, t) = ^ w n (r±, t) , 

n=l n=l 

s(r,t) = $> n (r ± ,t), (45) 



in which Nf is the number of filaments and r± = y/x 2 + y 2 . Note that the representation (45) 
is rather general and includes as well the case of no filamentatory structure, when Nf — 1. 
The number of filaments Nf will be defined later in a self-consistent way. From the point 
of view of quantum field theory, the appearance of spatial structures corresponds to the 
existence of nonuniform field vacuum [53,54]. In the case when different filaments are not 
correlated with each other, one has 

u m (r±,t)s n (r±,t) t)s n (r ± ,t). (46) 

This condition is typical of the turbulent regime, when the filaments are not mutually corre- 
lated [36-42] . Since filaments do not interact with each other, their location in the transverse 
cross-section is random. The radiation inside each filament is mainly concentrated along the 
filament axis, fading away at the distance much larger than the filament radius. In general, 
there can simultaneously exist the filaments of different radii. 

Let us consider an ra-th filament. And let the radiation in this filament be an order of 
magnitude weaker at the distance R n from its axis than at the latter, so that the intensity 
function w n (R n , t) is an order of magnitude smaller at r± = R n than w n (0, t) at r± = 0. The 
effective radius of the filament, r n , can be defined by the averaging relation 

TjT / w n (r ± ,t)r ± dr ± = w n (r n ,t) . (47) 

If the profile of the intensity function w n is of normal law, that is, 

w n (r ± , t) = w n (0, t) exp ^- , (48) 

where the filament radius r n plays the role of the standard deviation, then the relation (47) 
yields 

r n = j^j- 4 = 0.55 R n . (49) 
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All radiation of a filament, with an effective radius r n , is concentrated inside the envelop- 
ing cylinder of radius R n . Let us define the averaged functions 

u(t) = — ( u n (r±,t)dr, w(t) = ^- f w n (r±,t) dr , 

s(t) = — [ s n (r ± ,t)dr, (50) 

where the averaging is over the enveloping cylinder of volume V n = nR^L, and the enumer- 
ation index n for the left-hand side functions is omitted in order to simplify the notation. 
For what follows, we shall need the definition of the coupling functions 



sin(/c r — kz) 

and 



a(t)= loP [ Q(ct-r) ^y kz ) dr (51) 
Jv n k Q r 



m=l*( 8(c« - r) COB "f " kz) dr . (52) 
Introduce also the averaged stochastic field 

m = ^r U{v,t)e- ikz dv (53) 

and the field potential 

= -id- Exe-^ + eW- (54) 

Then, substituting the representation (45) into the evolution equations (40) and averaging 
according to Eq. (50), we come to the equations 

^ = -i(u + (5s)u - (72 - as)u + fi-s , = ~ 2 (72 - as)w + {u* }\ + fiu)s , 

at at 

ds 1 

_ = -aw - ^{u*fi + fiu)s - 7i(s - C) , (55) 

describing the dynamics of each of the filaments. 

In this way, the representation (45), together with the averaging procedure (50), have 
made it possible to pass from the equations (40) in partial derivatives to the ordinary dif- 
ferential equations (55). Following further the scale separation approach [5,47-50], we can 
more simplify Eqs. (55) by taking into account the existence of small parameters related to 
the standard situation, when the attenuation rates are essentially smaller than the transition 
frequency, 

— <1, — <1, — <1. (56) 

UJq UJq UJq 

This tells us that the function u in Eqs. (55) is fast, as compared to the slow functions w 
and s, which are the temporal quasi-invariants of motion. The collective width 

r = 72 - as , (57) 

collective frequency 

D, = uj + ps , (58) 
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and effective detuning 

5 = lu -n = A - f3s (59) 

are also slow functions in time, as compared to the fast function u. The latter can be found 
from the first of Eqs. (55), which gives 

V 6 + iFj 5 + iT 

+s C at'y-^+^'Ut 1 . (60) 

Jo 

Without the loss of generality, it is possible to choose the phase of the external field E x so 
that Uq6. ■ Ei be real, that is, 

u* Q vi = u ul , u = u(0) . (61) 

The solution (60) has to be substituted into the second and third of Eqs. (55), which are 
the equations for the slow functions. Then the right-hand sides of these equations are to be 
averaged over time and over the stochastic variables according to the rule 



1 r 

lim - / < . . . > dt , 



keeping fixed all temporal quasi-invariants. The stochastic variable £, describing local fluc- 
tuations, is assumed to be such that 

« £(t) » = . (62) 

Define the effective attenuation 

f = T» + J^(l-e- r <), (63) 

in which \5\ < \Y\ and the first term is caused by the quantum local fluctuations resulting in 

the quantum attenuation 

73 = Re lim i £ dt J* < C(t)Z(t') > e -(^+ r )(*-*') dt' . (64) 

Thus we obtain the equations for the guiding centers 

dw ~ ds ~ 

_ = -2(7 2 - as)w + 2Ts 2 , — = -aw - Fs - 7i(s - C) ■ (65) 

UjL (JjL 

Here the fast field fluctuations have been averaged out, and only the slow dynamics of 
filaments is left. 
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5 Transition Dipolar Waves 

An essential part of the quantum local fluctuations, leading to the quantum attenuation 
(64), is due to the dipolar vector potential (27) entering the local field (29). In order to 
better understand the physical origin of these fluctuations, let us consider the dipolar part 
of the local field (29) having the form 

2A; d • A dip (r 4 , t) = i £ [b^Sr(t) - Cij S+(t)] , (66) 
where we use the notation 



a/3 



2! ^ V 'J / 7 ' J 9-7T 

Z7F a/3 Z7F a/3 



/ G(ct - \Ti - r'|) Fjj exp (-iA; |r, - r'|) dr' . 

J 1*7. F 



Leaving in the pseudospin equations (16) only the terms related to the dipolar part of the 
vector potential, we have 

~^ = -\siY. {m - aSt) . (67) 

These equations describe the pseudospin fluctuations, which can be characterized by the 
deviations 

5S- = S-- <S~>, 5Sl = St- < St > (68) 

from the corresponding average values. Linearizing Eqs. (67) with respect to the deviations 
(68), under the condition < S~ >= 0, shows that 

S~ = 5S~ , St = const . 

Employing the Fourier transforms for the pseudospin operators 

S? = E S k exp (ik • r,) , S^ = ^Y1 S i ex P H k ' '*) ( 69 ) 
k Iy j 

and, similarly, for the coefficients 

K = E h k eX P ' r i) » °ij = E °k ex P ( ik • r i) > 

k k 

where = — rj, we get the equations for S^ and S£ which read 

= -ifikSk - i\kSt , ^ = in* k S£ + i\* k S^ , (70) 
14 



where 

fi k = u -b* k < St > , \ k = c k <S z t > . (71) 
The solution to Eqs. (70) has the form 

S k =u k e-^ t + v* k e^ t , (72) 

with the spectrum 



u k = \J |/ifc| 2 - |A fc | 2 . (73) 

This means that the dipolar part of the vector potential generates local fluctuations realized 
as a kind of transition dipolar waves. Such waves are analogous to the dipolar spin waves in 
magnets [55]. The spectrum (73) is always positive, since \b k \ 2 = \c k \ 2 <C uoq. This implies 
that the transition dipolar waves are stable. In the long- wave limit k — > 0, one has bo = Co, 
because of which 

lim uj k = uo . (74) 

The transition dipolar waves are responsible only for local fluctuations, but they do not 
participate in collective effects which are related to the coupling functions (51) and (52). 
These functions are zero at the initial time t = 0, though they grow very fast. Collective 
effects come to play after the interaction time r int = a/c, where a is the nearest-neighbour 
distance and c light velocity. For a ~ 1CT 8 cm and c ~ 10 10 cm/s, this time is very short, 
r int ~ 10 -18 s. After this time, the coupling function (51) quickly reaches its maximal value 
#7 2 and (52) grows to g ,r y 2 - Here the coupling parameters 

7o f sin(k r - kz) 
g = p— dr 75 

72 Jv n k r 

and, respectively, 

g> = p^ I C0S(A f ~ kz) dr (76) 
72 Jv„ k r 

are introduced. Recall that for each filament there are its own coupling parameters, whose 
values depend on the filament characteristics and the number of atoms in the filament. 
Strictly speaking, collective effects appear already for two atoms [56,57]. However noticeable 
radiation coherence develops only when a large number of atoms is involved. 
For t ^> r int , the collective width (57) becomes 

r = 72 (1 - gs) (77) 

and the collective frequency (58) is 

VI = uj + g' l2 s . (78) 

When the seed field is negligibly small, so that \v\\ 72, then the effective attenuation (63) 
simplifies to 

f ~ 73 . (79) 
Therefore, for times after Tj nt , the evolution equations (65) rearrange to the equations 

dw 2 ds 

— = -272(1 - gs)w + 273s , — = -#72^-73S-7i(s-C) , (80) 
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characterizing the filament dynamics. Recall that Eqs. (80) describe the slow dynamics, 
since the uncorrelated fast oscillations have been averaged out in the process of derivation of 
Eqs. (80). The random fast fluctuations are of less physical importance, being, in addition, 
much smaller than the correlated slow functions [36-42]. 



6 Temporal Dynamics of Filaments 

Equations (80), describing the filament dynamics, are to be complemented by the initial 
conditions w = w(0) and s = s(0). The first dynamic stage in the interval < t < r int 
can be ignored because of a very short interaction time Tint- In order that at the initial time 
t = the coherence intensity would be an increasing function of time, such that dw/dt > 0, 
it is necessary and sufficient that the inequality 

^(gso - l)w + 73S0 > (81) 

be valid, which is a criterion of laser generation. Since the second term in Eq. (81) is always 
non-negative, a sufficient condition for laser generation is gs > 1. 

A necessary condition for the formation of filaments is the appearance of coherence 
between radiating atoms. At the beginning, there exists a quantum stage of spontaneously 
radiating atoms, before coherence develops to a noticeable amount. The incoherent quantum 
stage lasts till the crossover time t c , after which coherent effects become dominant, and 
filaments are being formed. During the quantum stage < t < t c , the development of 
coherence is caused by the quantum fluctuations related to the term 273s 2 in the first of 
Eqs. (80). At this quantum stage, the system is yet uniform. Filaments are well formed 
after the crossover time t c . The formation of filaments inside a uniform atomic system, due 
to local quantum fluctuations, is somewhat analogous to the formation of galaxies from a 
uniform matter, due to local fluctuations [58], or to the stratification in quantum systems 
[59,60]. But the principal feature of laser filamentation is that the atomic matter in the laser 
cavity does not stratify as such. Nonuniformity happens for the distribution of radiating 
atoms, that is, it occurs on the level of photons, but not atoms themselves. This is why the 
filamentation in lasers can be called the photon filamentation. 

To estimate the crossover time t c , we need, first, to consider the quantum stage of evolu- 
tion, when there is no yet self-organized coherence. Assuming that at the initial time there 
is no coherence imposed by external fields, so that w = 0, we get the equations for the 
quantum stage, 

dw o ds , , 

— = 2 l3 s\ — = -( 7l + 73 ) s + 7l C. (82) 

The solution for the population difference is 

s = ( s -^r~] exp{-( 7 i + 73) t} + 



7i + 73 / 7i+73 
At short time, when (71 + 73 )£ <C 1, this yields 

s^s -(s - -^-) ( 7l + 73 )t + ±( ao - _Z^_) ( 7 i + 73 ) 2 t 2 . (83) 



7i + 73 



) (7i + 73) t + Us ^—) (7i + 7 3 ) 2 t 2 - 

/ 2 V 71 + 73/ 
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The latter form, since usually 71 < 73, is valid if 73^ « 1. At this quantum stage, the 
solution for the coherence intensity, given by the first of Eqs. (82), essentially depends on 
the initial value s of the population difference. For an arbitrary s . we have 

w ~ 2 73 s 2 t _ 27 3 s (7i + 73) (so -^-) t 2 + 

\ 71+73/ 

+| 73(7i + 7 3 ) 2 (so - -If-) (2*0 - -If-) t 3 • (84) 
3 V 7i + 73 / V 7i + 73 / 

When s 7^ 0, then the leading term in the coherence intensity is linear in time, 

w ~ 2 l3 s 2 t (s ± 0) . 
But if s = 0, then Eq. (84) gives the cubic dependence 

w ^ 2 ~ 7i 2 73C 2 * 3 (so = 0) . 

During the quantum stage, the evolution of coherence is mainly due to the term 273s 2 
in the first of Eqs. (80). But at the following coherent stage, the term 2 r y 2 (gs — l)w in this 
equation becomes dominant. Therefore, the crossover time t c can be defined as that one, 
where the two terms of different nature coincide, 

l2 {gs - l)w = 73s 2 (t = t c ) . (85) 

Note that under the initial condition s = 0, no noticeable coherence can evolve. This 
follows from Eq. (85) which, together with Eqs. (83) and (84), yields t c ~ T 2 = I/72. The 
rising of coherence should occur during the crossover time f c < T 2 . 

If s 7^ 0, then Eq. (85) gives the crossover time 

*c = — t 7, ^ • (86) 

71 (SO - C) + l2{gSQ - 1)Sq + 73S0 

Usually, one has 71 72 ~ 73, and, as is discussed above, one should have gs > 1 in order 
that collective effects could be important. Then the crossover time (86) can be reduced to 

To 

t c = — f— _ . (87) 

2(<?s -l) V ; 

For gso ^> 1, it is evident that t c <^ T 2 . At the crossover time (87), solutions (83) and (84) 
can be approximated as 

w(t c ) ~ 2 73 t c s 2 , s(t c ) ~ s . (88) 

After the crossover time t c , coherent effects become important. For the coupling param- 
eter g ^> 1, keeping in mind that 72 ~ 73, one has (772 ^> 73. In the standard situation, 
7i *C 72, because of which in the temporal interval t c < t T\ = I/71 equations (80) can 
be simplified to the evolution equations 

dw ds 

— = -272(1 - gs)w , — = -g^ 2 w , (89) 
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describing the transient coherent stage of filament dynamics. These equations can be solved 
exactly, resulting in the solution 



«;=PH sech 2 (^ 

\9l2 I T P , 



s = - — tanh f 

9 912 I t p , 



(90) 



whose integration parameters are obtained from the initial conditions (88). The pulse width 
7 P is defined by the relations 



7p = lg + 2(^72) 2 73^ c So , lg = (9 s o ~ l )l2 , l P r p = 1 
The pulse time, taking into account that 7 3 t c <C 1, reads as 

T 2 



Tn 



And the delay time is 



gs ~ 1 



t = t c + -| In 



73^ 2 gp 
G/*o - l) s 



Ip + lg 



(91) 



(92) 



(93) 



7p lg 

Recall that solution (90) is valid only if the pulse time r p is much less than T\. When r p is of 
order or much longer than T 1; as in many experiments [36-42], than one needs to deal with 
the total Eqs. (80). 

In view of the fact that 7 3 t c <C 1, the pulse width can be written as 



1p - (gs - 1)72 + 



9 2 l2l3t c sl 



gso-i 

where it is assumed that gs > 1. Then the delay time (93) becomes 

2(gs - 1) 2 



to = t c + t| In 



g 2 i3t c sl 

The latter, under strong coupling, when gs ^> 1 and r p ~ 2t c , rearranges to 

2 

to - t c + t c In 



(94) 



l2,tc 

with t c ~ T 2 /2^ So . 

At the final dynamic stage, when t >> 7\, we have to consider all terms in Eqs. (80). 
These can be written in the form 



~dt 



V2 , 



(95) 



with the notation 

vi = -272(1 - gs)w + 273s 2 , w 2 = -g^iw - 13s ~ii(s -() . 
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The stability properties of the solutions to Eq. (95) are characterized by the Jacobian matrix 
J(t) = [Jij(t)], with the elements 

Ju ^ ~dw = 2l2 ^ 9S ~ l > ' ~ds = l2gw + 73S ' 

J21 = = ~9l2 , J22 = -q^ = -7l - 73 • (96) 

Considering the stationary solutions to Eqs. (95), given by the equations v± — V2 — 0, 
we will present only those of them that correspond to stable fixed points. It is possible to 
separate out three cases, depending on the values of the coupling g and pumping parameter 


When g( <C — 1, then the stationary solutions 

w*^^f\, s* ~ C ( 1 - -^-) , (97) 

i2\g\ V 7iKI/ 

present a stable node, since the eigenvalues of the Jacobian matrix, that is, the characteristic 
exponents, are 

J\ - -7i - 7s , J2 ^ -272I5CI • 

Then the sole transient coherent pulse occurs, after which the radiation intensity diminishes 
to a value proportional to w*. 

For \g(\ <C 1, the stationary solutions 

2 



* / 7iC \ 73 
w ~ — 

\7i + 73/ 72 



l 7i(7i -7s)^C 



s * „ lA ( 1 _ ) m 

7i+73 V (7i + 73)V 



(7i + 7s) 2 

also corresponds to a stable node, with the characteristic exponents 

J\ ^ -7i - 73 , h ^ -272 • 

Hence again, only a sole transient pulse can appear. 

Finally, when both the coupling and pumping are strong, so that g( ^> 1, the stationary 
solutions 

w* ~ ^ , s* ~ - ( 1 - 1 (99) 

725 5 V TiffC/ 

represent a stable focus, with the characteristic exponents 

Jl,2 - - ^ (7l + 73) ± iWoo , 

where the effective asymptotic frequency 



^00 = \J2gClil2 ■ 

This means that the strong pumping supports the pulsing regime of radiation, when a set 
of bursts arise, being, at long times, separated from each other by the asymptotic period 



^00 v ac 
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The pulsing regime of radiation resembles pulsating stationary states happening in some 
statistical systems [61]. 

In this way, the flashing dynamics of each radiating filament depends on the related 
coupling parameter g and the level of stationary pumping defined by the pumping parameter 
(. In general, there are the following stages of evolution. The interaction stage < t < r int , 
when each atom radiates independently, before the signal reaches its nearest neighbours. 
This stage is very short and usually can be neglected. The quantum stage Ti nt < t < t c , 
when the radiation dynamics is mainly governed by local quantum fluctuations. After the 
crossover time t c , coherent collective effects become important, signifying the presence of 
the coherent stage. At this latter stage, depending on the value of the pumping parameter 
C, there occurs either a sole coherent pulse or a series of coherent bursts. If a stationary 
pumping is absent, then the coherent stage lasts in the interval t c < t < t + r p , and for t ^> t 
it changes to the relaxation stage, when the coherence intensity exponentially diminishes to 
the low level (97). 



7 Spatial Structure of Filaments 

Filaments are randomly distributed in the transverse cross-section of the laser cavity, evolving 
in space and time independently of each other. The characteristics of each filament essentially 
depend on the value of the related coupling parameter (75). For cylindric symmetry the latter 
can be presented in the form 



g = 2np — [ r ± alr ± [ 
72 Jo J- 



fL/2 sin(/c y / r5_ + z 2 — kz) 



-L/2 



dz . 



(101) 



Keeping in mind the resonance condition k ps k and introducing the variable x = k(\Jr\ + 
z), we have 



P7o f Rn , f kL sin a; 
g = 2tc - — / rj_ arj_ / ax . 

K7 2 JO Jkr 2 . IL X 



(102) 



Since, according to the inequality (43), we have A < L, then the upper limit kL in the 
integral (102) can be replaced by kL — > oo. This gives 

' kr 2 
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k^ 2 Jo 



— — bl 

2 



r ± dr L , 



(103) 



with the integral sine 



Si(x) 



x S1B.U 7T 

du = — h 

o u 2 



x smu 



oo U 



du . 



Introducing the notation 



Trig 
XL 



(104) 



varying in the interval < tp < ttF and playing the role of an effective Fresnel number for a 
given filament, we transform Eq. (103) to 



9(<P) = * 



PloL 
k 2 ^ 2 



7T(f 



2<p 



Si(x) dx 



(105) 
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g '^) = -n^ rCi(x)dx, (106) 



In the same manner, the coupling (76) can be reduced to 

/c 2 72 Jo 

with the integral cosine 

Ci(x) = / du . 



u 

Performing the integration in Eqs. (105) and (106), we may write 
9^) = 71 \- ni P ~ 2y?Si(2y?) + 1 - cos(2y?)] , 

Thus, the coupling parameters are functions of the effective variable (104), which, in turn, 
depends on the enveloping radius R n related to the effective filament radius r n by Eq. (49). 
In general, the filaments of different radii can arise. However, some of them are more stable 
than other, because of which the overwhelming majority of filaments possess the radii close 
to a typical value. 

The distribution of filaments with respect to their radii and, hence, the typical radius, 
can be found by invoking the general method of probabilistic pattern selection [18,21,22]. 
Following this approach, we define the probability distribution 

p ^ ,t) = W) exp { -x ^'*)> (107) 

for a filament characterized by the variable tp at the moment of time t. Here 

X((p,t) = Re /*Tr %>,£') dt' (108) 
Jo 

is the expansion exponent, expressed through the Jacobian matrix J of the evolution equa- 
tions, and 

Z(t) = J exp{-X(ip,t)} dip 
is the normalizing factor. The expansion exponent (108) defines the local expansion rate 

A(p,t) = ^X(p,t) . (109) 

The latter can be represented as the sum of the local Lyapunov exponents [18,21,22]. The 
partial sum of only positive Lyapunov exponents defines the entropy production rate [62,63], 
hence, the latter does not coincide with the local expansion rate (109). 

Thus, the probability for the appearance of filaments, characterized by the parameter 
ip, is given by the probability distribution (107). As is evident, the most probable is the 
filament with a typical <p satisfying the principle of minimal expansion [18,21,22] 

maxp(tp,t) mmX(ip, t) <^=^> mmA((p,t) . (HO) 
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This general principle follows from the minimization of the pattern information [22] and can 
be employed for arbitrary dynamical systems. 

The problem of turbulent photon filamentation is described by the evolution equations 
(80) or (95). The corresponding Jacobian matrix is given by Eqs. (96) from where 

Tr%>,*) = - 7l -7 3 - 2 72 (l-«? S ), 

with g = g(ip) and s = s(t). For t ^> 7\, the expansion rate (109) can be presented as 

A(<p, t) ~ -7i - 73 - 2 72 (1 - gs*) . (Ill) 

Using Eqs. (97) to (99), we get 

A(y?, t) ~ - 7l - 73 - 2 72 (1 + \g(\) (g( < -1) , 

A(V9,t)^- 7 i-73 - 27 2 fl- (Kl«l), 

V 7i + 73/ 

A^, *) ~ - 7l - 7s - ^ (K»l)- (H2) 

The stationary pumping parameter £ is in the interval — 1 < ( < 1, depending on the level 
of pumping. When there is no stationary pumping, ( = — 1. One says that the pumping 
is weak, if — 1 < ( < 0, and it is strong, if < ( < 1. Keeping in mind that the coupling 
parameter g is positive, we see that there exist two different cases, when the stationary 
pumping is weak or absent, ( < 0, and when it is strong, ( > 0. According to the principle 
of minimal expansion (110), the minimum of the expansion rate corresponds to the maximum 
of g(ip) if ( < 0, and to the minimum of g(^p), if £ > 0. The extrema of g(<p), as follows from 
Eq. (105), are given by the equation 

Si(2v?) = \ ■ (113) 

In the standard situation of absent or weak pumping, ( < 0, we have to look for the 
absolute maximum of g(tp). Then Eq. (113) gives ip = 0.96. From relation (104), we have 
R n = 0.55\/AL, and from Eq. (49), we find the typical filament radius 



r f = 0.3V XL . (114) 

The number of filaments can be estimated as Nf pa R 2 /R^, which yields 

N f pa 3.3F . (115) 

The linear dependence of the filament number on the Fresnel number is characteristic of the 
turbulent photon filamentation. 

Note that under strong stationary pumping (( > 0), when we need to look for the 
minimum of g((p), we would have ip = 2.45, hence, R n = 0.88^/\L and 77 = 0.5V AL. 

The formula (114) for the typical filaments can be compared with the radii observed in 
experiments. Thus, in different vapour lasers [25-29], one has 77 pa 0.01 cm. For the C0 2 
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laser and dye lasers, it was found [36-42] that 77 ~ 0.1 cm and 77 ~ 0.01 cm, respectively. 
All these data are in good agreement with formula (114). 

It is also worth mentioning that a similar kind of turbulent photon filamentation could 
arise in another type of matter, called photon band-gap materials. These materials possess 
a prohibited band gap, where light cannot propagate. The spontaneous radiation of atoms, 
with a frequency inside the prohibited band gap, is strongly suppressed [64]. However, if the 
density of doped atoms is sufficiently high, coherent interactions may develop (see review [5]). 
Then atoms can start radiating even inside the prohibited band gap. An unrealistic model, 
with the radiation length A ^> L much larger than the system size was considered [65] , where 
atoms with the resonance frequency at the band edge could produce collective spontaneous 
emission. Collective phenomena in the true atomic radiation of wavelength A <C L, with 
the atomic frequency inside the prohibited band gap were also considered and the effect of 
collective liberation of light was predicted [66-70]. Coherent radiation, accompanying this 
effect, should be realized by means of a bunch of turbulent filaments. 

In conclusion, the theory of turbulent photon filamentation in large-aperture lasers has 
been presented. This theory makes it possible to describe the spatial filamentary stricture 
as well as its dynamics. The typical filament radius, predicted by the theory, is in good 
agreement with experiments for different lasers. 
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